Symmetric excitation and de-excitation of a cavity QED system 
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We calculate the time evolution of a cavity-QED system subject to a time dependent sinusoidal 
drive. The drive is modulated by an envelope function with the shape of a pulse. The system 
consists of electrons embedded in a semiconductor nanostructure which is coupled to a single mode 
quantized electromagnetic field. The electron-electron as well as photon-electron interaction is 
treated exactly using "exact numerical diagonalization" and the time evolution is calculated by 
numerically solving the equation of motion for the system's density matrix. We find that the drive 
causes symmetric excitation and de-excitation where the system climbs up the Jaynes-Cummings 
ladder and descends back down symmetrically into its original state. This effect persists even in 
the ultra-strong coupling regime where the Jaynes-Cummings model is invalid. We investigate the 
robustness of this symmetric behavior with respect to the drive de-tuning and pulse duration. 
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I. INTRODUCTION 

A quantum two level system (TLS) interacting with 
a single mode of a quantized electromagnetic field is a 
central topic within the scope of circuit quantum electro- 
dynamics (QED). Typically, the light-matter interaction 
is weak enough to warrant the use of some version of 
the Jaynes-Cummings (JC) model [l[ to describe both 
static and time dependent systems [5-8]. However, 
recent progress in the field of circuit QED has enabled 
the fabrication of ultra-small mode- volume cavities where 
the light-matter interaction strength can be a consider- 
able fraction of a cavity photon energy. In this regime, 
evidence of the breakdown of the JC-model (with the 
rotating wave approximation) has been observed in su- 
perconducting [i| and semiconductor systems [ill HH ■ 

In previous work, we have gone beyond the TLS 
approximation and solved the many-body Schrodinger 
equation exactly for electrons embedded in a semicon- 
ductor nanostructure, subject to a single mode quan- 
tized EM field and an external classical magnetic field 
[T3 |. The electron-electron and photon-electron interac- 
tions were treated exactly using exact numerical diago- 
nalization (for details see Ref. [HI]). We predicted the 
failure of the JC-model (including anti-resonance terms) 
at high coupling strengths for a static and closed system. 
Here, we expand on that work by adding an explicit time 
dependence to the total Hamiltonian and investigate its 
dynamical properties. 

We investigate the effects of time dependent addition 
to the total Hamiltonian which does not depend on the 



cavity photon creation or annihilation operators. In the 
language of the JC-model, this means that we are per- 
turbing the atomic term of the total Hamiltonian as op- 
posed to the cavity field or interaction term. The time 
dependent term is a sinusoidal drive which is modulated 
by an envelope function which varies slowly compared 
with other characteristic time scales of the system. 

The paper is organized as follows. In Sec. |H] we de- 
scribe the static part of the Hamiltonian which contains 
the geometry of the semiconductor nanostructure as well 
as the electron-electron and photon-electron interaction. 
In Sec. IIIII we add a time dependent drive to the total 
Hamiltonian and introduce the time dependent observ- 
ables we are interested in and how we calculate them. 
Results and concluding remarks are presented in Sees. 
IIVI and [V] respectively 



II. DESCRIPTION OF THE STATIC SYSTEM 

We split the static part of the system's Hamiltonian 
rio into four terms, 



rlo — ri e + "Hem + Hp + H a 



(1) 
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They are the electronic part of the Hamiltonian (H e ), the 
cavity field Hamiltonian (Hem), the paramagnetic (H p ) 
and diamagnetic {rid) electron-photon interaction terms. 
The electronic part can be written as 

H e = H w + Hcoul (2) 
= J2 E i4di + \ X^il^Coui|rs)44cW r , (3) 

i ijrs 

where H w is the Hamiltonian of a finite quasi-one- 
dimensional (QID) quantum wire with hard walls at x = 
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Figure 1. (Color online) Charge densities for the SES's |1) (a) 
and 1 2} (b) in arbitrary units for zero magnetic field. Dark 
is low density and bright is high, (c), sinusoidal pulse with a 
Gaussian envelope function with a = 10 which gives 10 com- 
plete oscillations in the time interval r. (d), spatial variation 
of the drive potential r(a;) for /3 X L X = 3 (dashed red line), 
/3 X L X = 4 (thin green line) and f3 x L x — 5 (thick blue line). V 
is homogeneous in the y-direction so a 3D plot is not needed. 



bracket) with energy E^. The states \fi) and |/i) are 
connected via a unitary transformation |/i) = V|/i). We 
calculate V by diagonalizing H e in the basis of the eigen- 
states oiH w . Where needed, we use the notation |/i)jy to 
denote the /z-th electronic state containing N electrons. 
For example, |4)2 is the fourth lowest two electron state. 
Note that we have a closed system so the number of elec- 
trons is a conserved quantity. 

The cavity field EM Hamiltonian can be written as 
Hem = fudpCt'a where (r and a are bosonic creation and 
annihilation operators of a cavity photon with energy 
hiOp. The photon energy is typically chosen such that 
it matches a transition energy in the quantum wire (i.e. 
the photon field is near a resonance). 

By taking the long wavelength approximation for the 
quantized EM field, its vector potential can be written 
as Aem = ^EMe(a + a') with e a unit vector pointing 
in the field's direction of polarization. The paramagnetic 
interaction term can then be written as 



Hv — £< 



E 



d]d j g lJ (a + a r ) 



(5) 



where we have defined the electron-photon coupling 
strength £ c = qA^M^wO-w, which is the characteristic en- 
ergy scale for the photon-electron interaction. We have 
also introduced an effective dimensionless coupling tensor 
(DCT) 



9ij 



(6) 



±L x /2 and a parabolic confinement in the y-direction 
with characteristic energy MIq, containing several non- 
interacting electrons. The wire is subject to an exter- 
nal classical magnetic field B = V x A oxt = Bz. With 
the magnetic field, the characteri stic ener gy in the y- 
direction is modified to Ml w = U^VLq + us%, where lu c is 
the cyclotron frequency w c = qB/m*, with q the positive 
elementary charge and m* the electron effective mass. 
A natural length scale is a w = \J h/(m*£l w ). d\ and di 
are fermionic creation and annihilation operators of the 
single electron eigenstates (SES) \i), with energies Ei, 
which include the effect of the external magnetic field. 
Figs. QJa) and[TJb) show the charge density of the two 
lowest SES for zero magnetic field. 

"Hcoui contains the effect of the Coulomb interaction 
with the kernel 



y/(x - x') 2 + (y - y') 2 + V 2 



(4) 



where r\ is a small convergence parameter to regulate 
the singularity at r = r'. We denote the many electron 
eigenstates of 1L W as (we use Latin indices for single- 
electron states and Greek ones for the many-electron 
states) with the energy E^. These states are simply 
Slater determinants in the SES's We denote eigen- 
states of T-L e as (denoted now with a rounded right 



defining the coupling of individual single-electron states 
\i) and | j) by the photonic mode where ir = p + qA. ext is 
the mechanical momentum. It is useful to generalize g^j 
to 



Qixv = ^9ij{lA4 d i\ v ) 



(7) 



so that Gfj,v is the DCT which defines the coupling of 
many-electron states |/k) and \v) by the photon field. 
From Eq. we see that Q^ u depends on the geometry of 
the quantum wire and on the magnetic field. Note that 
the action of di and d\ is only known in the {|//)} basis 
(of non-interacting Fock states). To calculate for exam- 
ple di\/j.) we have to use di\(j,) = diV\fx) . This definition of 
Q^ v makes comparison with the JC-model easy, since the 
coupling energy £ jc in the JC model (in the £ jc (T x{a+<l < ) 
term) is related to the DCT via £jc = |£7 k a|£c where 
| re) and |A) are the two states chosen for the TLS ap- 
proximation (active states) [d|. We will refer to £ c as 
the electron-photon coupling strength (or simply the cou- 
pling strength) and £jq as the effective coupling strength. 
The diamagnetic interaction term can be written as 



£ 2 



— (aW + aa] 



(8) 



where Af e is the electron number operator. 
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Now that we have defined all the terms in ((T|), we can 
put the results together and expand Ho in the basis, 



U =y^E /J ,\fj,)(/j 1 \ + hupa^ a + ^2S c g^\fi)(v\(a + a f ) 

[A [IIS 




where N e is the number of electrons in the system. 

Next we proceed to expand Ho in a basis containing a 
preset number of photons {|M)} (eigenstates of a' a with 
eigenvalue M). We have then expanded Ho in the com- 
plete orthonormal basis {\a}} = {\p) <2> |Af)}. From the 
diagonalization process we obtain a unitary transforma- 
tion hi which satisfies lA\a) = \a) where \a) are eigen- 
states of Ho satisfying Ho\a) = E a \a). The states |a) 
are then used as a basis in time dependent calculations 
which are covered in section Mil 

III. TIME DEPENDENT HAMILTONIAN 

Now we add a time dependent drive to the Hamiltonian 
in the form (in first quantization) 

W(t) = WoT{x, y) sm(u> t)F(t) , (10) 

where Wo is the drive amplitude which has units of en- 
ergy, T(x, y) contains the spatial dependence of W(t) 
and F(t) is an envelope function for the sinusoidal term 
sin(a;oi). We choose u>o such that hajo matches some tran- 
sition energy in the quantum wire (i.e. the drive is near 
a resonance). We will refer to hu>o (u> ) as the drive en- 
ergy (frequency). In this work we will choose F(t) to be 
a pulse which varies slowly on the timescale 2tt/ojo and 
satisfies F(±oo) = and F(t) < 1. An example of such 
an envelope function is a Gaussian 

F(t) = e' it - to)2/T2 , (11) 

where r 3> 2tt /loq- It is useful to define the quantity a = 
ujot/(2it) which gives the number of oscillations within 
the time interval r (see Fig. He)). Note that a larger 
value of a translates into a longer pulse (not a faster 
oscillating pulse). As for the geometric part we will use 

T(x, y) = - e -0l^-^/2? + e -0U*+L./2)* ; (12) 

so T(x,y) = T(x) is a sum of two Gaussians of opposite 
sign that are centered on the opposite end of the quan- 
tum wire (see Fig.[]Td)). The reasons for this choice of T 
is that it couples the time dependent part of the Hamil- 
tonian strongly to a cavity field that is polarized in the 
x-direction. 

To add W(t) to Ho, we need to calculate the second 
quantization generalization of W(t) in the {\p)} basis 



using 

w{t)=^m\m4dAv)\ti{v\ 

ijfj,u 

=W wa.(u t)F(t)^2T IJU ,\fi)(v\ (13) 

where the matrix element T ^ characterizes the coupling 
of individual Coulomb interacting many-electron states 
|^i) and \v) by W{t). In ([T3]), W(t) is simply a unit op- 
erator in the photon Fock space since it does not depend 
on a or at 

Now that we have the matrix representation of the to- 
tal Hamiltonian T~L = Ho + VV(£), we can calculate the 
time evolution of the system by integrating the equation 
of motion 

ih^ = [H,p}. (14) 

where p is the density matrix of the system. This is done 
numerically using a Crank-Nicolson method. We let the 
system start out in the ground state and investigate its 
excitation by W(t). 

The observables we are interested in are the mean num- 
ber of photons (N)(t) and energy E{t) = (H)(t) which 
can be calculated using (A)(t) = Tr [p(t)A(t)} where A 
is some observable. 



IV. RESULTS 

Throughout this section, we use L x = 300 nm, MIq = 
1.0 meV, f3 x L x = 4, a = 90, m* = 0.067m e and 
e = 12.4eo (GaAs parameters). Unless otherwise stated, 
both the cavity photon energy hui p as well as the drive en- 
ergy huio are on resonance between the electronic states 
|1) and |2) with de-tuning S = 10~ 4 Ai2, where A12 = 
E-2 — Ei is the energy difference between the two low- 
est electronic eigenstates. We let the system start out in 
its ground state and always use an envelope function of 
the form in (fTTj) . Where needed, we will differentiate the 
de-tuning of the cavity field and the drive energy as 5 P 
and Sd respectively. We will only consider x-polarization 
for two reasons. First, it couples the states |1) and |2) 
strongly (large |^i2|). Second, the quantum wire is ap- 
proximately (for small magnetic field) harmonic in the 
^-direction. We want to compare our results to simpler 
TLS models and for a TLS approximation to be appli- 
cable, there needs to be some anharmonicity present to 
minimize excitation to states outside of the TLS Hilbert 
space. 

Fig. [5] shows energy spectra for Ho with one and two 
electrons. From the figure we can see 4-6 of the lowest 
polariton pairs in the JC-ladder, along with other states 
that are not a part of the TLS approximation. Note how 
much denser the two electron spectrum is. The range 
of coupling strength is similar to what we use for time 
dependent calculations later on. 
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Figure 2. (Color online) (a), Energy spectra for one electron 
at zero external magnetic field as a functions of the effec- 
tive coupling strength £jc = \Gi2\£c where \Gi2\ = 0.300 and 
Aia = 0.197 meV. (b), same as in (a), but for two electrons 
with |(?i2 1 = 0.676 and A i2 = 0.521 meV. Note that the en- 
ergy spectrum has been shifted so that the ground state has 
zero energy at zero coupling strength. 



Fig. [3] shows the excitation energy E cxc (t) = E(t) — 
E(t = 0) of the system as a function of time for var- 
ious driving field amplitudes (Wo)- The most interest- 
ing aspect of the figure is the symmetric excitation and 
de-excitation of the system. What is happening is that 
the system is climbing up the Jaynes-Cummings ladder 
and descending back down symmetrically. This would 
not be surprising for a low coupling strength where both 
the anti-resonant and diamagnetic terms of the electron- 
photon interaction Hamiltonian can be ignored. In that 
case, when both the drive and photon energy are on reso- 
nance, it is possible to find quasi stationary states which 
are periodic in time (see Ref. @). The effect of the enve- 
lope function is then to adiabatically tune the system's 
quasi ground state. The adiabatic nature of the envelope 
function prevents transitions beyond this slowly shifting 
quasi ground state. When the envelope function goes 
to zero again the system is adiabatically shifted to its 
original state. Like mentioned earlier, this behavior is 
expected at low coupling strength, however, for the plots 
in Fig. [31 the electron-photon coupling is rather large 
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Figure 3. (Color online) Excitation energy vs time for a sys- 
tem containing one electron at zero magnetic field. The en- 
velope function F(t) is a Gaussian of the form in (|11[1 with 
to = 2r and a = 90, giving 90 oscillations in the time inter- 
val r. The effective coupling strength is £jc/(hup) — 0.13. 
The system starts out in the ground state at t = 0. The 
different curves are for different Wo (see Eq. (|13p ~) where 
W0/A12 = n x 0.428. 



(£jc/(7iWp) — 0.13). In fact, for this coupling strength, a 
visible deviation from the JC-model is apparent in Fig. 
[5Ja) where the polariton splittings don't show the char- 
acteristic linear ±y/n£jc splittings which we expect from 
the JC-model (for further comparison with the JC-model, 
see Ref. [l2| or [H]])- Another interesting aspect of Fig. [3] 
is that the value of the maximum energy -E max at t ~ 2r 
scales quadratically with Wo . This will be addressed later 
in this work. 

Fig. 2] shows the excitation energy of the system as a 
function of time for various effective coupling strengths 
(£jc), keeping the drive amplitude (Wo) constant at 
W0/A12 ~ 3.4. Note that in £ JC = \Qi 2 \£c, only £ c 
can be varied since depends only on geometry. The 
same excitation and de-excitation pattern is apparent, 
even for the ultrastrong coupling £jc/(Sw p ) = 0.512. We 
did not perform calculations for higher coupling because 
of issues with numerical convergence at higher coupling 
strengths (see Ref. [13] for detailed convergence calcula- 
tions). Counter- intuitively, the maximum excitation en- 
ergies at t ~ 2r in Fig. [5] is proportional to the inverse 
square of the coupling strength £ c . We can therefore see 
that the maximum excitation energy follows the relation- 
ship £ max oc (W /£ c ) 2 - 

By doing calculations for many different values of Wo , 
£ c and magnetic field strength, we were able to deduce 
the relationship 



E,, 



A(N e 



\r 12 \w 



(15) 



\Qi2\£c 

Where A(N e ) is a dimensionless constant that depends 
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Figure 4. (Color online) Excitation energy vs time for a sys- 
tem containing one electron at zero magnetic field. All param- 
eters are the same as in Fig. [3] except £jc is varied instead of 
Wo- Note that changing the electron-photon coupling shifts 
the ground state energy as can be seen in Fig. [2] That is 
why the above energy is shifted by the ground state energy 
Eo. The three different numerical values 0.128, 0.256 and 
0.512 are the value of Ejc/i^p) used for for the correspond- 
ing curve. 

only on the number of electrons but not on the geometry 
of the system. Deducing the (Wq/£ c ) 2 relationship was 
straightforward since Wq and £ c are both parameters we 
can freely change. However, to show the |r 12 /(?i2| 2 rela- 
tionship we had to vary the matrix elements Ti2 and Q12. 
We did this by doing calculations for different magnetic 
field strengths (see Fig. [5]). Varying the magnetic field 
changes the geometry of the system and thereby changes 
both T\2 and Q12. By doing this we eliminate the possi- 
bility of the symmetric excitation de-excitation behavior 
being a special case of the geometry we have chosen or 
due to the time reversal symmetry (which the magnetic 
field breaks). Note that changing the magnetic field and 
electron number alters the energy spectrum so A12 is 
not always the same. From Fig. [5] we see that the agree- 
ment with Eq. (|15p is very good for one and two electrons 
while large deviation is apparent for three electrons. The 
reason for this is that for more electrons, the energy spec- 
trum is much denser and previously inactive states start 
to have bigger influence. Least square fit with (|T5|) gives 
A{\) = 0.26, A{2) = 0.21 and A(3) = 0.058 which shows 
that a larger electron number translates into a weaker 
response to the drive. 

We find that the mean number of photons follows the 
same relationship as the energy 

where A^ max is the peak value of (N)(t). The constant of 
proportionality A(N e ) is the same as for (fT5|) . 
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Figure 5. (Color online) Maximum excitation energy for many 
values of £ c , Wo and B with £ c and Wo in the range of that 
of Figs. [3] and [4] and B in the range 0.0-0.5 T. Results are 
shown for one (a), two (b) and three (c) electrons with 52, 
47 and 18 data points respectively (some data-points overlap 
and are not distinguishable on the graphs). The blue lines 
are least square fit with |15]) giving A(l) = 0.26, A{2) = 0.21 
and ^4(3) = 0.058. The results are very robust for one and 
two electrons but start to show significant deviations for three 
electrons. 

Because the geometric dependence of Eqs. (fT5")) and 
P6p is encoded in the two matrix elements Ti2 and Q12 
it is tempting to say that a TLS description is appli- 
cable and we can ignore the inactive electronic states 
(all electronic states except |1) and |2)). This is true 
for small coupling strength where the JC-model can be 
used and the diamagnetic interaction term is small. How- 
ever, as can be seen from Fig. [6j a TLS description 
starts to fail at coupling strengths of £jc/(fiw P ) ~ 0.1. 
What is surprising in Fig. [5] is that for coupling strength 
£jc/(fiw P ) ~ 0.1, including the diamagnetic interaction 
term in a TLS model gives much less accurate results 
than if it is ignored. This behavior is observed when 
£jc/(fu*Jp) > 0.10. It seems that the effects of the dia- 
magnetic interaction term are canceling the effects of 
the inactive states. However, for even higher coupling 
strengths, the TLS model (with or without the diamag- 
netic interaction term) will fail because it can't account 
for the complicated effect of inactive states in an energy 
spectrum such as the ones in Fig. [5] 

We will end this section by investigating the robustness 
of the results. Fig. [7] shows results with different values of 
the drive de-tuning. We find that symmetric excitation 
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Figure 6. (Color online) Excitation energy vs time. All 
three plots are for one electron, 5 = 0, Wo = 0.24 meV, 
£ c = 0.08 meV and \G 12 \ = 0.3, giving £ JC /(Hu p ) ~ 0.13. 
Curve 1 shows the exact results using our full model. Curve 
2 shows results obtained using a TLS approximation using 
only two states for the electron Hilbert space basis where the 
diamagnetic interaction term is ignored. Curve 3 shows re- 
sults using a TLS approximation including the diamagnetic 
interaction term. Note how the TLS approximation with dia- 
magnetic term drastically overestimates the excitation energy 
and how the system is not symmetrically de-excited. 
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Figure 7. (Color online) (a), excitation energy vs time for two 
values of drive de-tuning for one electron with £jc/(fot) P ) — 
0.13, Wo/£ c = 5 and B — 0. The two curves are for a de- 
tuning <5d/Ai2 = 0.02 (red-dashed) and <5d/Ai2 = 0.05 (solid 
blue), giving A ~ 5.0 and A = 1.0 respectively. The symmetric 
behavior is lost in the latter case, (b), Same as in (a) except 
Bic/iJujJp) — 0.26 and Wo/£ c = 3. The two curves are for 
a de-tuning <5<j/Ai2 = 0.01 (green-dashed) and <5d/Ai2 = 0.1 
(solid black), giving A ~ 14 and A ~ 1.4 respectively. The 
symmetric behavior is lost in the latter case. Note that in 
both (a) and (b), the cavity photon de-tuning is held constant 



means that we need at least 10 complete oscillations 
within the time period t. This is the same value of a 
which is used for illustrative purposes in Fig. QJc). 

Throughout this work we used a Gaussian shaped en- 
velope function. It is worth noting that the above men- 
tioned symmetric behaviour is not a special case for this 
choice. We tested this by varying the envelope function 
such that F(t) = cxp { — \t — to\ n /T n }, with n in the range 
1-5. This change does affect the shape of the excitation 
energy curves such as those in Figs. [3] and 2] but the ex- 
citation and de-excition is still symmetric and has the 
same maximum value as with a Gaussian shaped enve- 
lope function. 

Finally, we note that so far we have only used envelope 
functions that are symmetric around t — to- We found 
that relaxing this restriction does change the symmetric 
exciation and de-excitation behavior such as not being 
symmetric about t ~ to anymore, but as long as the 
envelope function is slowly varying, the system still ends 
up in the same state as it started in (the ground state). 
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Figure 8. (Color online) Excitation energy vs time for differ- 
ent values of a. (a), one electron, B = 0.2 T, £ c — 0.16 meV, 
Wo = 0.48 meV, A i2 = 0.172 meV and f JC /(fiw P ) = 0.245. 
(b), two electrons, B — 0.1 T, £ c = 0.08 meV, Wo = 0.32 meV, 
A12 = 0.511 meV and £ 3C /(huj p ) = 0.101. In both (a) and 
(b), we see that in order observe the symmetric behavior, we 
need to have a > 10. Note that the time axis has been shifted 
in both (a) and (b) such that the pulse maximum is at t — 0. 



CONCLUDING REMARKS 



de-excitation behavior persists as long as 

We also investigated what is the minimum pulse length 
needed to observe the symmetric behavior (see Fig. HJ). 
We find that the minimum value a ~ 10 is needed. This 



We have described a method to compute the time evo- 
lution of a many-body, multi level, Coulomb interacting 
electronic system which is coupled to a single-mode quan- 
tized EM field. The Hamiltonian contains explicit time 
dependence on the form of a sinusoidal drive which is 
modulated by an envelope function. The envelope func- 
tion is a Gaussian shaped pulse that varies slowly com- 
pared to characteristic time scales of the system. Both 
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the drive frequency and the cavity photon frequency are 
on resonance with a transition in the electronic system. 
The only approximations which are made are the finite 
size of single- and many-body bases as well as the finite 
precision of the numerical integration. The convergence 
with respect to these parameters is carefully controlled. 

We observed a symmetric excitation and de-excitation 
of the system where the drive causes the system to climb 
up the Jaynes-Cummings ladder. When the pulse dies 
out, the system descends down the ladder back to its 
original (ground) state. This is not a simple adiabatic 
effect because even though the envelope function varies 
slowly, the drive frequency matches that of the cavity 
photons as well as a transition in the electronic nanos- 
tructure. 

We found that when the system exhibits symmetric 
excitation and de-excitation, the maximum expectation 
values of the energy and photon number follow a simple 
relationship which is independent of the specific geome- 
try, but does depend on the number of electrons in the 
system. The effect is very robust for one electron, slightly 
less so for two electrons and considerable deviation is 
observed for three electrons. This is due to the much 
denser many electron energy spectra, where states which 
are not a part of the Jaynes-Cummings ladder start to 
have greater influence. 

We compared the results of our full model to that of a 
TLS model. We did this by doing calculations where only 
two states are used as a basis for the electronic part of 
the Hamiltonian. We found that for the relatively weak 
coupling strength £jc — 0.13, including the diamagnetic 



interaction term gives much worse results than if it is 
ignored. The results of the TLS model without the dia- 
magnetic term is in agreement with the results of the full 
model. 

We tested robustness of the above mentioned sym- 
metric behavior with respect to the magnitude of the 
drive de-tuning. We found that the maximum allowed 
value of the relative de-tuning was 1-2%. Note that 
this value depends strongly on geometry, coupling- and 
drive amplitude (see Eq. Q17p). Finally, we investigated 
how long the pulse needs to be to see the symmetric 
effect. We found that approximately 10 complete 
oscillations are needed within the time period r, where 
r characterizes the duration of the pulse via Eq. ([TTj) . 
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